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Abstract 

In this paper the Randall- Sundrum model with brane-localized curvature terms 
is considered. Within some range of parameters a compact extra dimension in this 
model can be astronomically large. In this case the model predicts small deviation 
I/-) . from Newton's law at astronomical scales, caused by the massive modes. The existence 

of this deviation can result in a slight affection on the planetary motion trajectories. 
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1 Introduction 

m ■ 

Nowadays models with brane-localized curvature terms are widely discussed in the literature. 
In paper £Q it was argued that matter on the brane can induce a brane-localized curvature 
\q • term via quantum corrections, which appears in the low-energy effective action. An attractive 
feature of this model is a modification of gravity at ultra-large scales, which can be very 
! interesting from the cosmological point of view. Unfortunately a strong coupling effect was 
found in the model Nevertheless, the DGP-proposal can be utilized not only in the 

models with infinite extra dimension. For example, in EJ Ej some models with compact 
extra dimension and brane-localized curvature terms were discussed. The model proposed 
in [H] possesses some very interesting feature - an extra dimension in this model can be 
Ph! astronomically large. In the present paper we discuss another model with astronomically 
large extra dimension. It is based on the Randall-Sundrum solution for the metric jB] and 
was discussed in [7] in the case of sub-millimeter extra dimension. Here we consider this 
model from another point of view - in the case of a very large size of extra dimension. 
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2 The setup 

The action of the model considered in |7] has the following form 

S = S g + Si + S 2 + S in d, (1) 

where S g , Si, S 2 and S in d are given by 

S g = -^J (R-A)j=g#xdy, (2) 

51 = -j^J(R-A)6Mj=§*xdy, 

52 = -^^J{R-A)5(y-R)^d 4 xdy, 
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fiind is a dimensionless parameter and A = —12k 2 . Here g^ v is the induced metric on the 
branes and the subscripts 1 and 2 label the branes. The model possesses usual Z 2 orbifold 
symmetry. We also note that the signature of the metric guN is chosen to be (— , +, +, +, +). 
Obviously, the model admits the Randall-Sundrum solution for the metric 



ds< 



^MNdx* 4 dx N = / y fll/ dx tJ 'dx u + dy l 



(3) 



where 7^ = e 2r7 ^r]^ l ,, rj^ u is the Minkowski metric and the function a{y) = —k\y\ in the 
interval — R < y < R. The parameter k is positive and has the dimension of mass. The 
function a has the properties 



d^a = —k sign(y), 



dla 



-2k(S(y) - 8(y - R)) 



-2k5. 



(4) 



We denote k = V 167rG, where G is the five- dimensional gravitational constant, and 
parameterize the metric jmjv as 



gMN = IMN + kh 



MNi 



(5) 



Hmn being the metric fluctuations. In [7] the linearized equations of motion (for the field 
Iimn) corresponding to the action (JTJ) were derived. It was shown, that in the gauge 



h 



0, h. 



hAA(x) 
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these linearized equations possess an additional symmetry under the transformations 

1 / 1' 



h^(x,y) 



d^d v ip{x) 



(6) 



(7) 



<j>(x)^Kj>{x) + ip(x), (8) 

which do not belong to the gauge transformations. The validity of the gauge (JOJ) was thor- 
oughly checked in 00113 • It is evident that with the help of transformations (JTJ), (JBJ) we can 
impose the condition 4>{x) = 0. In other words it is necessary to examine the next order of 
perturbation theory to get the equation for the radion field. But since the radion field does 
not interact directly with matter on the brane, it seems that the linear approximation is not 
destroyed by the radion. The key feature of the model is that the only physically relevant 
case, in which tachyons and ghost are absent (and the symmetry (J7J), (jHJ) is preserved) is 
when the matter (and the induced term S in d) exists on brane 1 only. Brane 2 (at the point 
y = R) can be interpreted as a "naked" brane, i.e. a brane without matter on it (see [Zj). 

The gravity on the brane in the linear approximation is described by the field h^, which 
satisfies the following effective equation of motion [Zj 



e~ 2 °U 



2k 



h 



[IV 



+ d'ih^ ) - 2k 2 V + 2kh^5 



(9) 
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t^v denoting the energy-momentum tensor of matter on the brane 1 and t = fft^. Using 
this equation, we can get the equation for the zero mode of h^, which has the form h Q = 
Let us multiply equation Q by e 2<J and integrate it over coordinate y. Using the 



2a 



orthonormality conditions, which have the form [7j 



dye 



-2a 



\Kv) + l*(y -R) + %^%) 



In ■ 



(10) 



we get 



□ a 



kk 



(n; 



A fully analogous procedure was made in the case of RSI model in jHj. 

In the next section we will estimate the effects produced by the massive modes. We will 
not solve equation ©, as it was made in for the RSI model, we will thoroughly estimate 
the masses of the modes, their wave functions and coupling constants to matter on the brane. 
We would like to remind that we will use other values of parameters, not that used in [7j. 



3 Wave functions of the massive modes 

It was shown in |7| , that equation for the wave functions of the massive modes has the form 

I (e- 2(J D V + dl V) " 2A; 2 V + 2 ^V - ^'^^ + ( 12 ) 

+ ^%)e- 2ff D V = 0. 

Following the footsteps of [7||H], we arrive at the relations: 

* n (y) = N n (N (t n e kR ) J 2 (t n e^) - J {t n e kR ) N 2 (t n e~ CT )) , (13) 

where N n is a normalization constant, and 

N (t n e kR ) Jo (*„) - Jo {t n e kR ) N (t n ) + (14) 
+ n 2 md [No (t n e kR ) J 2 (t n ) - Jo (t n e kR ) N 2 (t n )] = 0. 

Here h^ u {x,y) = J2^ n (y)h^ u (x), ^ n (y) is the wave function of the corresponding massive 
mode, t n — ^ and m n is such that Oh™(x) = ra^/i" (x). Equation (JT4*jl defines the mass 
spectrum of the theory. 

Let us choose fli n d such that fli n d » 1 (we will specify its value later) and kR « 1. For 
the relatively large values of ? one can use the expansions for the Bessel functions for large 
values of arguments (see, for example, [HQ]). In this approximation masses of the modes are 
defined by 

(e^-l)!^-™, (15) 
where n is integer. Small corrections to these values are such that 

*" (0) ^W (16) 
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This formula can be obtained from equation (|14j) (see Appendix 1). It seems that this 
form of the wave functions (inverse proportional to m^) can be obtained only if there exist 
brane-localized curvature terms. It is necessary to note that since the main purpose of 
[7j was to examine the structure of linearized equations of motion, only a brief analysis of 
^"(0), resulted in a very rough estimate \l/ n (0) ~ was carried out in that paper. This 
estimate is sufficient for not rejecting the model in the case of sub-millimeter extra dimension. 
Equation shows that contributions of the massive modes appear to be suppressed much 
stronger, and for the case of sub- millimeter extra dimension (an example in j7j) an effective 
theory on the brane actually describes massless 4-D gravity. At the same time the form of 
(fTT)|) is very important for the case of very large extra dimension, which will be considered 
below. 

It is not difficult to show that integral in (|1U|) is always positive and N n ~ ^| for relatively 
large m n (see Appendix 2). Thus the form of interaction of the massive modes with matter 
on the brane looks like 

\ [ d'x (k"£* n (0K„W + tf^t- E -K»W u ) > ( 17 ) 

2 Jbrane \ ^ MpAnd £j 71 J 

where a n ~ 1, M Pi = VL ind / \/ Gk and j is such that for n > j one can use expansion 
formulas for the Bessel functions with good accuracy. One can see, that coupling constants 
of the massive modes are suppressed by the factor VL in d and strongly depend on the number 
of a mode. It is not difficult to show that for n < j the factor k^/ n (0) is of the order of 
1 / (QindMpi) . This estimate can be obtained from (flUjl. (|T3|l and (fT4^) using the fact that the 
Bessel and Neumann functions are of the order of unity or lesser for the masses, which are 
the solutions of (|14|) : and the fact that combination in the square brackets in (fT4*j) is equal 
to \l/ n (0) (accurate to the factor N n ). 

Now let us discuss the possible values of parameters of the model. We would like to 
mention, that effective theory exactly of this type appears in the model considered in [3] , so 
we will use some results obtained there. Because of the fact that coupling of Kaluza-Klein 
modes to matter on the brane is defined by the factors (\l/ n (0)) 2 and the sum ^ \ is rapidly 
convergent, one may not to worry about the strong affection of massive modes because of 
their quantity (it was thoroughly checked in [5]). This shows an importance of the fact 
that (pH)|) is inverse proportional to and iV n is proportional to m n . Let us consider R 
of the order, say, of the size of the Sun. For this case constraints on Qmd following from 
the astrophysical experiments are very severe (see 0) - fl in d > 10 5 (in our notations). The 
multi-dimensional Planck mass M should be of the order of 10 GeV to get the proper four- 
dimensional gravitational constant, which is defined by (jllj) . Thus the hierarchy problem is 
solved in an analogous way to that proposed in fU^]. I n this case the massive modes can 
become apparent on the distances of the order of extra dimension's size. It can lead to some 
slight affection on the Mercury motion. 

Of course, the size of extra dimension can be much lesser: for example, from the size of 
the Earth down to the sub-millimeter scales and even smaller. In this case constraints on 
the value of Qmd following from the astrophysical data can be substantially relaxed, and the 
only claims are that the values of G and VL in( i with given R should give the proper value 
of four- dimensional gravitational constant and correspond to the constraints obtained from 
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the top-table gravitational experiments. Because of suppression of the massive modes, the 
phenomenology of this model differs considerably from that one of any standard "large extra 
dimensions" scenario (this issue was discussed in detail in [3]). 

4 Discussion and final remarks 

In this paper we constructed a self-consistent model with astronomically large compact 
extra dimension. One may ask: why we use warped geometry, whereas the same effect can 
be achieved in the model with flat background 0? Firstly, in the paper [5 3 the radion 
field (44-component of the metric fluctuations) was taken into account not throughout all of 
calculations. It is well-known that existence of this field can change the theory considerably, 
as it was in the case of DGP-model [21 13] • It is not probable that the same problem can arise 
in the model described in j^j, nevertheless it is necessary to solve corresponding equations 
of motions with regard to the radion field more accurately, checking the validity of imposed 
gauge conditions. In our case considered above the radion field is absent in the linear 
approximation and one can forget about it (and about the problems with ghosts, see (21 
UniED at least at the classical level. Secondly, the main advantage of warped geometry is 
that it ensures the proper tensor structure of the massless graviton propagator (in the case 
of flat geometry this property seems to be forbidden by the equation for the 44-component 
of metric fluctuations). It is also necessary to note, that the value of fl inc i was chosen to 
correspond an appropriate value in More realistic minimal value, corresponding to the 
experimental data, can be a little bit larger (because of the warped background, see (JEJ)). 
But it is evident that this difference should not be considerable. Some other interesting 
features of this setup were mentioned in 

Now let us discuss our choice for the size of extra dimension. Solution of equations of 
motion for the /ii/-component of metric fluctuations includes some term (see [7]), which is a 
pure gauge from the four-dimensional point of view of an observer on the brane. This term 
is proportional to the factor p-, and in the case of very large size of extra dimension (and 
correspondingly very small k) the linear approximation brakes down. Exactly the same thing 
happens in massive gravity at the small graviton mass limit, see for example ^2]- Thus the 
size of extra dimension of the order of the Sun's size seems to be "safe" from this point of 
view. For a cases of much larger size of extra dimension it is necessary to carry out a more 
detailed analysis including examination of the non-linear effects. 

We would like to mention, that only the case e kR ~ 1 was considered in this paper. 
The case e kR » 1 seems to be quite interesting, because there can appear an additional 
suppression of contributions of the massive modes by the exponential factor, which could 
relax the constraints on Qi n d- At the same time the lowest modes will have lower masses 
~ ke~ kR in this case. Anyway this situation calls for more detailed investigation. Also it 
would be interesting whether the symmetry analogous to (JTj), (jH)) could arise in other models, 
for example, with infinite extra dimension, which admit modification of gravity at ultra-large 
scales. 
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5 Appendix 1 

Let us calculate estimates (fTHjl and (fT7)j) for the case ? >> 1 and kR « 1. For large values 
of arguments we can use the following expansions for the Bessel functions [10 



N (t) 




id 
id 



cos ( t — 



71 



1 / 7T 

— sin [t 

St V 4 




. 7T\ 1 / 7T 
sin [t COS \t 

4/ 8t V 4 





7T \ 15 / 7T 
COS t 7T Sin t 7T 

4 / 8t V 4 



. , 7T \ 15 / 7T 
Sin t 77 H COS t 71 

1 4 / St V 
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(19) 
(20) 
(21) 



Substituting JIED-JSU) into JUJ), we get 

sin ((A - l)t) + 
= -nl d (sin((A-l)t)- 

where t = ^, X = e kR . For the relatively large t solutions of this equation in the zero 
approximation are defined by 

U 



1 1 
8t ~~ 8\1 _ 
15 1 
si + 8X1 



Tin 



cos ((A - l)t) 
cos ((A - l)t) 



A - 1 



where n is integer. Corrections At n to t n at the first order are defined by the equation 

(-l) n At n (A - 1) (Q 2 ind + l) = (-1)' 
from which one gets 



_ J_ + J_ + + is^L 



3At n 8At n 



8t 



At 



(A - l)t n 

Finally for the function \l/ n (0) we have 



15 1 
— + 



16 



# n (0) « JV r , 



vrt n V Trt n X V 8f2? nd t n 



16 



•I)" 



4JV„ 



-i) 
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6 Appendix 2 

Let us show, that 

dye 



-2a 



i - \Kv) + \Kv - r ) + % W 



W 2 > o, 



if $ satisfies the equation 



(22) 



- {e- 2a m 2 m + <9 2 *) - 2k 2 ^ + - 
2k 2k yy > 



(23) 



Here 



^2(1/) 



^2(1/), 

Mt 6 



fefl 



/ / m - 



Jo( ¥ e 



„T ( m 

N <T e 



N is normalization constant. It is easy to show that 



R 



dye-^Zi(y) 



R 



m 

1 = T e 



2k f tR 



tZ%(t)dt, 



(24) 



where t = f, 



tR = ir e • With the help of the Bessel equation for the function Z 2 (t) 



Z 2 (t) + 



z 2 (t) 



t 2 



+ 1 - 3 = 



^2(0 = 4fZ 2 {t), one can get for the integral 



2A; 



to 



2fc /"*« 
2 



to 



ti 2 2 (t) + 



^(t) 

4Z 2 (t) 



AZ 2 (t) 



Z 2 (t)-tZ 2 (t) 



dt 



t 



dt + 



2k 



-tZ 2 (t)Z 2 (t) 



tR 



t 



where integration by parts were made. Using the boundary conditions for Z 2 (t), which follow 
from the equation (}2"3*|) and have the form 



tZ 2 (t) + 2Z 2 (t)--Z 2 (t) 



t=t R 



t 2 



i 1 



tz 2 (t) + 2z 2 (t) - -z 2 (t) + n 2 nd -z 2 (t) 



0. 



0. 



t=t 



we can show that integral (}2*2*j) is equal to (accurate to the factor N 2 ) 

Az 2 (ty 



2k ( f tR 

2 



m 



tZ 2 {t) + 



t 
7 



dt+2 [Z 2 (t)] 



(25) 



Substituting the recurrent formula for the Bessel function 

Z 2 (t) = Z 1 (t) - jZ 2 (t), 

where 

Z 1 (t) = N (^e kR ) J, (t) - Jo (je kR ) N, (t) 

into (J2SJ), we get 



( f R \tZl{t) - AZ 2 {t)Z 2 {t)] dt + 2 [Z%(t)] \*A = — 2 f* tZ\{t)dt. (26) 
\Jt 1 J J m j to 



2k 

m 2 

It is evident that this integral is always positive. Thus the integral in (J2"2"|) is also positive 
for any real N. 

Now let us calculate (J2T?)) for the relatively large m. Using the formulas analogous to 
(fTS|) - (f2Tj) (see [TO] ) one can easily show that 

Z x {t)a l_ cos(t-t R ). 

Thus 



m 2 J to m 2 TT 2 tu m 2 7r 2 



and N ~ ^ for fci? w 1. 
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